Measuring the top Yukawa coupling at the ILC at y^s = 500 GeV 
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We report on the feasibility of the direct measurement of the top Yukawa coupling gt at the 
International Linear Collider (ILC) during its first phase of operation with a center-of-mass energy 
of 500 GeV. The signal and background models incorporate the non-relativistic QCD corrections 
which enhance the production cross section near the ti threshold. The e + e~ — > tiH signal is 
reconstructed in the 6-jet + lepton and the 8-jet modes. The results from the two channels are 
combined. The background processes considered are e + e~ — >• ibW~ /tbW + (which includes e + e~ — > 
ti), e + e~ — > iiZ, and e + e~ — > tig* — > ttbb. We use a realistic fast Monte-Carlo detector simulation. 
Signal events are selected using event shape variables, through jet clustering, and by identifying 
heavy flavor jets. Assuming a Higgs mass of 120 GeV, polarized electron and positron beams with 
(P e - , P e + ) = (—0.8, +0.3) , and an integrated luminosity of 1 ab -1 , we estimate that the e + e~ — > tiH 
events can be seen with a statistical significance of 5.2 a, corresponding to the relative top Yukawa 
coupling measurement accuracy of \Agt/gt \ = 10%. 

PACS numbers: 13.66.Jn, 14.65. Ha, 14.80.Bn, 
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I. INTRODUCTION 

The standard model (SM) of elementary particle 
physics stands on two pillars. The first pillar is the 
gauge principle, which has been verified by precision elec- 
troweak measurements. The second pillar consists of the 
electroweak symmetry breaking which is yet to be tested 
by experiment. The discovery of the Higgs boson will be 
of particular importance in explaining the mass genera- 
tion mechanism. Within the SM, the Yukawa interaction 
of the top quark and the Higgs boson generates the mass 
term which breaks the electroweak gauge symmetry The 
measurement of the strength of the top Yukawa coupling 
gt can shed light on the mechanism behind the generation 
of the top quark mass. 

The top-quark Yukawa interaction could be measured 
indirectly using the production mechanism of the Higgs 
boson through the top quark loop at the LHC experi- 
ments. The indirect measurement unfortunately cannot 
give a full description of the top-quark Yukawa interac- 
tion, for suppose we observe an anomaly in the Higgs 
production cross section; it would be difficult to distin- 
guish whether this effect is due to an anomaly in the top 
Yukawa interaction itself, or because there are contri- 
butions from unknown particles propagating in the loop 
connecting the initial state and the Higgs boson. In order 
to distinguish these two effects, it would be highly desir- 
able to measure the top Yukawa interaction directly. At 
the LHC, the direct production process gg — > tiH in the 



H — > bb channel is marred by jet combinatorial back- 
ground While the H — > 77 or H — > t + t~ channels 
are expected to yield cleaner signals [2[ , which could allow 
for the discovery of the gg — > ttH process, the uncertainty 
in the top Yukawa coupling value would be affected by 
the potentially large uncertainties in the Higgs branch- 
ing fraction measurements. We show that a future e + e~ 
linear collider, such as the International Linear Collider 
(ILC) , can play a critical role in the determination of the 
top Yukawa coupling through the direct measurement of 
e + e~ — > ttH in the H —> bb channel. 

Feasibility studies of the top quark Yukawa interaction 
at a future e + e~ linear collider have a long history 0,0]. 
A serious feasibility study of a direct measurement of 
the top Yukawa coupling using the process e + e _ — > ttH 
at the center-of-mass (CM) energy yfs = 800 GeV was 
performed in [5J, |6| , which incorporated realistic exper- 
imental conditions expected at a linear collider experi- 
ment. More recently, there have been increased interests 
in how well the top Yukawa coupling can be measured 
in the first phase of a linear collider experiment, whose 
CM energy reaches up to yfs = 500 GeV. An analysis for 
y/s = 500 GeV was carried out in 0. It was noted that 
at this energy the bound-state effects between t and t en- 
hance the ttH production cross section significantly, since 
the relative momentum of t and t is typically small [JhLJ] ■ 
A reanalysis was performed in the Snowmass workshop, 
incorporating the enhancement effect by ti resonance for- 
mation as well as an enhancement effect that can be ob- 
tained by polarizing the e + e~ beams [l5[. The conclu- 
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sion was that the top Yukawa coupling can be measured 
to roughly 10% accuracy, including statistical errors only, 
with an integrated luminosity of 1 ab _1 . 

In this paper, we investigate the feasibility of measur- 
ing the top Yukawa coupling at ^fs = 500 GeV using 
the process e + e _ — > itH for the Higgs mass of 120 GeV. 
The new aspects of this study as compared to the previ- 
ous ones are as follows. We implement the enhancement 
factor by it bound-state effects into the event generator, 
both for the itH signal and the itZ background events; 
the latter is particularly important since the expected 
measurement accuracy of the top Yukawa coupling is sig- 
nificantly affected by the number of these background 
events. In addition, we perform a fairly detailed detector 
simulation which takes into account the realistic energy 
resolution of the calorimeter components. (See Sec. IIV[) . 

This paper is organized as follows. In Sec. UH we 
present our method for including the enhancement by 
it bound-state effects. In Sec. IIII1 the signatures of the 
itH process and the possible background processes are 
outlined. The analysis framework used for the event 
generation and the detector simulations is discussed in 
Sec. IIVI We discuss the event selection procedure in de- 
tail in Sec. [V] for the 6-jet plus lepton mode analysis, 
and in Sec. IVII for the 8- jet mode. We summarize the 
accuracy estimate of the top Yukawa coupling measure- 
ment in Sec. I VIII The measurement is assumed to be 
dominated by the statistical uncertainty. 



II. INCLUSION OF tt BOUND-STATE EFFECTS 



Theoretical analyses of the it bound-state effects on the 



cross section for e + e 



ttH with the next-to-lcading 



logarithmic accuracy are given in [14| . To our knowl- 
edge, there have been no analyses of it bound-state ef- 



fects for the background process e 



ttZ which are 



included consistently with the e + e~ — > itH signal. In 
our event generator, all the QCD corrections are incorpo- 
rated consistently only with the leading-order accuracy. 
Hence, for consistency, we incorporate the it bound-state 
effects on the signal and background cross sections with 
the leading-order accuracy; we also incorporate some of 
the important next-to-leading order corrections in the 
bound-state effects. 

There arc a number of (tree-level) Fcynman diagrams 
contributing to each of the processes e + e~ — > itH, 



ttZ and e 



+ ,:- 



ttg*, where the t(t) subse- 



quently decays into bW + (bW~). Let us denote these 
amplitudes for the process i —> f as A t j(i —> /). The 
tree level amplitudes are modified as follows: 

A#(i -> /) = [A t - t {i -> /)] tree 

x V K i^f x F (s t t,p;m t ,r t ,a s ). (1) 

F represents a process-independent enhancement factor 
that incorporates it S'-wave bound-state effects; de- 
notes the CM energy of t and t as reconstructed from the 



final bW + bW~ system; p is the three-momentum of t in 
the CM frame of t and t; m t and r t denote the pole mass 
and width of the top quark, respectively. Close to the 
threshold of it pair production, this factor F incorporates 
the bound-state effects according to the non-relativistic 
bound-state theory, while for higher values of 
factor F is smoothly interpolated to unity: 



F 



G(E,p) 
1, 




2m t <C m t , 



(2) 



The non-relativistic Green function is defined by 



(E + iT t )-\ + V QC v(r) 

m t 



G(E,p) 



d 3 xe- l P- s G(E,x), 



(4) 



where r = \x\ and Vq cd (r) is the next-to- leading or- 
der QCD potential [14[. Go(E,p) is the non-relativistic 
Green function of a free it pair, which is defined via 
Eqs. ([3]) and ((4]) after setting Vqcd(?*) to zero. The 
enhancement factor F is explained in more detail in 
Ref. [H]. 

In Eq. ([T]) Ki-+f denotes a process-dependent hard- 
vertex correction factor, which arises as a part of the 
next-to-leading order corrections. To a good approxima- 
tion this factor is independent of kinematical variables 
for the signal process e + e _ — > itH . In fact with a choice 
Ki-yf = 0.843, we reproduce the e + e~ — > itH differen- 
tial cross section at s/s = 500 GeV shown in Fig. 5(a) 
of Ref. [l3|. (The next-to- leading-logarithmic curve in 
the threshold region of it and the 0(a s ) curve at higher 
or lower Eh-) We adopt this value of K^f for the 
signal process. For the background processes, we choose 
K^f = 1, since these factors are unknown and since 
the deviation of these factors from unity is part of the 
next-to-lcading order corrections not fully accounted for 
in our analysis. 



III. SIGNAL AND BACKGROUND PROCESSES 

The diagrams for e + e~ — > itH with subsequent top de- 
cays t — > bW + (t — > bW~) are shown in Fig. [TJ The first 
and second diagrams containing the top Yukawa coupling 
g t are the targets of this study. The contribution to the 
cross section coming from the third diagram, where the 
Higgs radiates off of the intermediate Z boson, is negli- 
gible at y/s = 500 GeV, as shown in Fig. [2] As a result, 
the number of itH events is proportional to gf to a very 
good approximation, which enables us to perform a sim- 
ple analysis by event counting. 

In this study, the Higgs boson is reconstructed in the 
two 6-jet mode H — >• bb (68%). Our itH signal can be 
classified into three groups, depending on the decay mode 
of the W bosons. Their branching fractions are 
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FIG. 1. Feynman diagrams for the e+e — > process. 

(i) 8-jet mode: 45%, 

(ii) 6-jet + lepton mode (e or fi): 29%, 

(hi) 4-jct + 2-lepton mode (ee, e/i, or /i/i): 5%, 

where we have omitted the contribution of the top decays 
to tau (t —> br+iy T and t —> br~V T ), since we only recon- 
struct electrons and muons from the top in this study. 
The 8-jet mode and the 6-jet + lepton mode are chosen 
for reconstruction. 

The following processes are identified as possible back- 
ground sources which can mimic the ttH signatures: 

(i) e+e- -> tbW-/tbW+ -> bW+bW~, 

(ii) e+e- -> tiZ -> bW+bW-bb, 

(iii) e+e- -> tig* -> bW+bW-bb. 

The cross sections for these processes are shown as a 
function of yfs in Fig. [5J We will refer to the e+e- — > 
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FIG. 2. Production cross section of the e+e~ — > ttH sig- 
nal (shown with and without tt bound-state effects) , together 
with those of the main background processes, ttH (Higgs radi- 
ated off the Z boson), tiZ, tt, ibW~ fibW+ (denoted as tbW), 
and tig* — > tibb, as a function of the CM energy without beam 
polarizations. The initial state radiation and beamstrahlung 
effects are included. 



state consists of up to two b jets, as opposed to four b jets 
for our ttH signal. The tbW channel can be therefore re- 
duced to a small fraction by identifying the flavor of the 
b quarks in the final state (^-tagging) and by counting the 
number of b jets. Because of the large tbW cross section, 
a significant amount of tbW background remains even if 
there is a small rate of event mis-reconstruction, which 
occurs equally likely for events in and away from the top 
pair resonance, thus making it important to include the 
non-resonant contributions. 

In contrast to the tbW process, the processes tiZ and 
tig* can have identical final states as those of the ttH 
process if the Z boson or the hard gluon g* decays into 
a bb pair. In this case, the signal extraction will de- 
pend strongly on the resolution of the Higgs mass re- 
constructed from the two &-jets. The unpolarizcd cross 
section for tiZ is 1.3 fb, including the ii bound-state 
effects similar to that expected for the signal process; 
without including this correction, the cross section be- 
comes 0.7 fb. For tig* — > tibb, the unpolarizcd cross 
section is 0.7 fb. We note that there is no ii bound- 
state enhancement in the tig* process because the ii 
system is not a color singlet in this case. The cross 
sections at y/s = 500 GcV for our signal and back- 
ground processes are summarized in Tab. |TJ The signal 



TABLE I. Cross sections at y/s = 500 GeV for the signal and 
background processes are shown for the different beam polar- 
izations. The e+e- — > ttH and e+e~ — > tiZ processes include 
the ti bound-state effects. The ttH, tiZ, and tig* processes 
all decay as bW+bW~bb while the ibW+ /ibW~ process (de- 
noted as tbW) decays as bW+bW~ . The number of events 
N used in this study is shown for each sample, along with its 
equivalent luminosity £.. 



Process 




a (fb) 


N 


£ (ah" 1 ) 


e L e R ~+ 


ttH 


1.07 


5.00 x 10 4 


47.8 


e z4 -> 


tiZ 


4.04 


5.00 x 10 4 


12.4 


e L e R ~+ 


tig* 


1.93 


5.00 x 10 4 


25.9 


e L e R ~+ 


tbW 


1633 


1.00 x 10 7 


6.1 


e R e l -> 


ttH 


0.45 


5.00 x 10 4 


92.6 


e R e L ~+ 


tiZ 


1.32 


5.00 x 10 4 


37.8 


e R e l -> 


tig* 


0.86 


5.00 x 10 4 


58.2 


e R et -> 


tbW 


700 


1.00 x 10 7 


14.3 



ibW~/ibW+ process as e+e" -> tbW. The e+t 



tbW 



and background samples have been produced with pure 
beam polarizations. Unless otherwise noted, our results 
weight these samples to match the beam polarizations of 
(P e -,P e+ ) = (-0.8,+0.3) [13. 



IV. ANALYSIS FRAMEWORK 

Signal and background events are generated using the 
physsim [l8[ event generator, based on the full helic- 
ity amplitudes including gauge boson decays, calculated 
using HELAS and BASES [2C)j], which properly takes 
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into account the angular distributions of the decay prod- 
ucts. For the event generation, the following values are 
used: a(M z ) = 1/128, sin 2 6 W = 0.230, a s = 0.120, 
M w = 80.0 GeV, M z = 91.18 GeV, M t = 175 GeV, and 
Mjj = 120 GeV. The effects of initial state radiation and 
beamstrahlung are included. The it bound-state effects 
results in a roughly twofold increase in the itH signal 
cross section at y/s = 500 GeV, as shown in Fig. [3] The 



TABLE II. Detector parameters, p.pr arid E are measured 
in units of GeV. The angle 9 is measured from the beam axis. 

Detector Resolution Coverage 

Vertex detector a h = 7.0 © (20.0/p sin 3/2 6) Jim \ cos6\ < 0.90 
Drift chamber ap T /P T = 1.1 x 10 _4 p T 0.1% | cos6>| < 0.95 
ECAL <J E /E= 15%/x/E © 1% | cos 6\ < 0.90 

HCAL a E /E = A0%/x/E®2% cos 6>| < 0.90 
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340 345 350 355 360 365 370 375 380 

m ft (GeV) 

FIG. 3. Differential production cross section as a function 
of the invariant mass of the it system for unpolarized beams 
(P e - ,P e +) = (0.0, 0.0) with y/s = 500 GeV. The shaded his- 
togram represents the leading-order values. The white his- 
togram includes the tt bound-state effects. 



four-momenta of the final-state quarks and leptons are 
passed as input to Pythia 6.4 [21j for parton showering 
and hadronization. The detector response is simulated 
using the QuickSim [22j fast Monte-Carlo detector simu- 
lator. 

The detector consists of the beam pipe, a vertex de- 
tector, a drift chamber, an electromagnetic calorimeter 
(ECAL), and a hadronic calorimeter (HCAL). The cross- 
ing angle of the beams are also taken into account. Each 
hit in the tracking detector is smeared according to the 
detector resolution specified in Tab.im For each charged 
particle, the parameters describing its helical trajectory 
are smeared according to the full covariance matrix of 
the parameters. Calorimeters arc simulated down to the 
level of individual cells with possible overlaps of energy 
deposits from nearby particles. Each hit in the calorime- 
ter cell is smeared according to Tab. [TTJ The calorimeter 
cell hits are clustered and then matched to the tracks of 
charged particles. ECAL clusters which are consistent 
with a charged track are subtracted based on the parti- 
cle flow approach (23|. HCAL clusters whose energy is 
consistent with charged hadrons are removed, while for 
clusters with inconsistent energy matching, as is the case 
when neutral hadrons are present, HCAL energy deposits 
are statistically subtracted by an amount weighted by the 
geometrical overlap between the charged hadrons and the 
HCAL clusters. 



V. ANALYSIS OF THE 6-JET + LEPTON 
MODE 

We describe the event selection for the analysis of the 
6-jet + lepton mode first, followed by the 8-jet mode. 
Similar techniques are used in both modes. The main 
differences between the two analyses are the presence of 
a lepton in the 6-jet + lepton mode and the number of 
jets in the final state. 



A. Identification of an isolated lepton 

In the 6-jet + lepton mode, the lepton from the W —> 
iv tends to be energetic and isolated from the rest of the 
event. To identify such a lepton (e or fx), a cone with a 
half-opening angle 9 cone is constructed around each track 
(lepton candidate). The cone energy E conc is defined 
to be the sum of the energy of all the tracks inside the 
cone, excluding the lepton candidate. The value of 9 CO nc 
which gives cos 9 conc = 0.98 is found to be optimal for our 
event selection. Fig. @] shows the distribution of the cone 




J Isolated leptons 
| Leptons in Jets 




60 80 100 

Lepton Energy[GeV] 

FIG. 4. Distribution of the cone energy and the lepton energy, 
shown for leptons from W decays (black dots) and leptons 
originating in heavy flavor jets (gray dots). The black curve 
shows the cut boundary for the lepton selection; leptons below 
the curve are identified as isolated leptons. 

energy versus the lepton candidate energy. The energetic 
isolated leptons from the W decay have a high lepton 
energy and a low cone energy, thus populating the lower 
right region, shown as black dots in Fig. 3J Leptons from 
heavy flavor jets are likely to be less energetic and have 
a higher cone energy, shown as gray dots in Fig. 01 The 
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selection of isolated leptons is performed by applying a 
cut on the cone energy which varies as the lepton energy 
and is given by the equation E conc < y/6(E( — 15) where 
£J C one an d Eg are given in units of GeV. For the 6-jet 
+ lepton event selection, we require one and only one 
isolated lepton. 



B. Event shape 

To exploit the differences in the event topology be- 
tween e + e _ — > ttH and e + e~ — > tbW events, we use the 
thrust T, defined as [H HI 



T 



Ei \ n -Pi 
max ' 

i™i =1 Li Pi 



(5) 



where the index i runs over each reconstructed parti- 
cle with pi corresponding to its momentum, and T is to 
be maximized with respect to the unit vector h corre- 
sponding to the axis in which the overall event shape 
is stretched. The thrust T tends to unity for 2-jet-like 
events while it tends to 1/2 for isotropic events. Because 
the it events tend to have fewer jets, with a higher av- 
erage jet energy, T tends to be higher for tbW events 
compared to ttH events. The requirement of T < 0.77 is 
found to be optimal in the 6-jet + lepton analysis. 



C. Jet clustering 

We employ the Durham jet clustering algorithm [26j 
to separate the event into 6 jets, after taking out the 
isolated lepton. In the Durham algorithm, each particle 
is regarded as a jet on its own to begin with; a jet pair i 
and j gets combined if the pair has the lowest Y%j value 
which is defined as 

,E]}(1 



Y — 



min{_E, 



cos 6; 



E 2 



(0) 



where 9ij is the angle between the momentum vectors 
of the two jets, and E cm = y/s. The jet clustering is 
allowed to continue until there are two jets remaining, 
with the value of being recorded at each transition 
from n jets to n— 1 jets, which we call Y n ^ n _i. It is found 
that Y5_).4 is useful for discriminating ttH 6-jct + lepton 
events from tbW events. This is due to the fact that, after 
identifying the isolated lepton I from the semileptonic 
decay t — > bW — > blv, tbW events cannot have more 
than four jets without a gluon emission. Taking all the 
backgrounds into account, we require Y5_>4 > 0.005. The 
jet configuration for n = 6 is used for the rest of the 
analysis of the 6-jct + lepton mode. 



D. Identification of heavy flavor jets 

The identification of jet flavor is critical for the sup- 
pression of tbW background due to the differences in the 



number of 6-jets in the final state between tbW events, 
which produce two 6-jets, and our ttH signal, which re- 
sults in four 6-jets. The jet flavor is identified by looking 
at the number of secondary tracks belonging to the jet. 
We count the number of tracks with an impact parameter 
significance (in three dimensions) greater than a certain 
threshold value Q; if the count is equal to or greater than 
a certain number Nq, the jet is identified as a 6-jet. The 
two numbers (Q,Nq) are optimized in our 6-tagging se- 
lection. 

In the 6-jet + lepton analysis, we define two criteria for 
the identification of 6-jets. We define loose 6-jets as jets 
passing the 6-tagging requirement of (Q,Nq) = (2.0,2); 
tight 6-jets are defined by those passing the 6-tagging 
requirement of (Q,Nq) = (2.5,4). Note that the set of 
tight 6-jets is a subset of loose 6-jets. We require at least 
four loose 6-jets in the event; two out of the four are also 
required to pass the tight 6-jet criteria. 

The 6-tagging efficiency is estimated using a sample of 
Z —> qq events at y/s = 91.2 GeV and is found to be 81% 
(loose) and 47% (tight). The rate of incorrectly identify- 
ing a jet originating from a lighter quark as a fake 6-jet is 
estimated to be 40% (loose) and 3.2% (tight) in a sample 
of c-jets, and 0.5% (loose) and 0.08% (tight) for s, d, and 
w-jet samples. In a multi-jet environment, the 6-tagging 
efficiencies decrease due to overlapping jets. For 6-jet 
+ lepton and 8-jet events, this effect typically reduces 
the 6-tagging efficiencies by roughly 10%. In principle, 
given the fact that the 6 — > c mis-tagging rate is rather 
large, final states involving charm quark jets such as the 
process e + e~ — > ttg* — > bW + t>W~cc can be a source of 
background. This effect is not included in our analysis. 
However, we expect that the use of a more sophisticated 
6-tagging method will reduce the charm contamination 
in the selected 6-jets (27j . 



E. Top and Higgs reconstruction 

The Higgs candidate {H — > 66) is formed by requir- 
ing one tight 6-jet and one loose 6-jet. The hadronic top 
candidate is formed by combining three jets, one of which 
must be (at least) a loose 6-jet. Because there are multi- 
ple possible ways to combine the six jets in this way, we 
define the quantity % 2 as 



x 2 = 



- M 



1-1 



OH 



Mi 



<7 W 



(7) 

and choose the jet combination which minimizes the x 2 
value. Here, rrijj {rrijjj) is the invariant mass of the 
two-jet (three-jet) system; the two jets ji and ji are 
used to form the Higgs candidate, while j'3, j'4, and j'5 
are the three jets used to reconstruct the top candidate 
which decays hadronically. The masses M t , My/, and 
Mfj are taken to be the same values used in the event 
generation. The widths a t , ay/, and o~h correspond to 
the mass resolutions in the case of perfect jet clustering 
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and jet combinations. These values arc determined to 
be a t = 14.3 GeV, a w = 9.3 GeV, and a H = 17.7 GeV 
by combining the reconstructed four-momenta of final 
particles in the e + e~ — > UH sample using Monte-Carlo 
information. 

Final cuts are applied to the resulting invariant mass 
distributions. For the 6-jet + lepton mode, we require the 
top mass to be in the range of 140 < m t < 205 GeV and 
the Higgs mass to be in the range of 95 < mjj < 150 GeV, 
where the range has been optimized in steps of 5 GeV. 



pendent from each other, allowing for a straightforward 
combination of the two results. 



B. Event shape 

Similarly to the 6-jet + lepton analysis, the thrust vari- 
able T is used to reduce the tt background. It is found 
that T < 0.7 is optimal for the 8-jet mode. 



F. Results 

We summarize the yields after applying each cut for 
the case of polarized beams (P e - , P e +) = (—0.8, +0.3) in 
Tab. IIII1 where the yields are normalized assuming an 
integrated luminosity of 1 ab . The resulting distribu- 
tions for the thrust, ^5^4, the top candidate mass, and 
the Higgs candidate mass, after applying all the other 
cuts, are shown in Fig. [5] The signal significance is esti- 
mated to be 3.7, corresponding to the measurement ac- 
curacy of the top Yukawa coupling of \Ag t /gt\ — 14%. 
With unpolarized beams (P e -,P e +) = (0.0,0.0), the sig- 
nificance becomes 2.9, corresponding to \Ag t /g t \ = 17%. 



VI. ANALYSIS OF THE 8-JET MODE 

A. Isolated lepton rejection 

In the 8-jet analysis, the 6-jet + lepton mode can be- 
come a source of background as a result of splitting the 
jets by the jet clustering procedure. To reduce this kind 
of background, we look for isolated leptons using the 
same prescription used in the 6-jet + lepton analysis. 
Events containing one or more isolated leptons are dis- 
carded in the 8-jct analysis. This procedure ensures that 
the 6-jet + lepton and 8-jet samples are statistically inde- 



C. Jet clustering 

In the 8-jet mode analysis, the jet clustering is per- 
formed over all particles in the event to form eight jets. 
We keep the jet transition values Y n -y n -i. The value for 
Ys->7 is found to be useful in discriminating ttH events 
from it events. We require Ys->7 > 0.0009 in the event 
selection. 



D. Identification of heavy flavor jets 

We follow a similar procedure as in the 6-jet + lepton 
mode for the identification of 6-jets. We use a different 
optimization for the tight 6-jet, which is modified to be 
(Q,N Q ) = (3.0,2). The definition of the loose 6-jet re- 
mains the same. 



E. Top and Higgs reconstruction 

The Higgs candidate (H — > 66) is formed by requiring 
one tight 6-jet and one loose 6-jet. One of the top candi- 
dates is required to contain a tight 6-jet, while the other 
top is required to have (at least) a loose 6-jet. Because 
there are multiple possible combinations of jets, we de- 
fine the quantity x 2 similarly to the 6-jet + lepton mode 
as 



x 2 = ( m hh -M H \ 2 + f m j3jij5 -Mt \ 2 + / m jejTh - MA 2 + f m hU -M w \ 2 + / m jejT - M w \ 2 ^ 



and choose the combination of jets which minimizes the 
X 2 value. Here, j\ and j'2 are used to form the Higgs 
candidate. The three jets j'3, j'4, and .75 are used to re- 
construct the first top candidate, while j6j 37 1 and j$ are 
used to reconstruct the second top candidate. The same 
values for at, aw, and an are used as in the 6-jet + 
lepton analysis. 

Final cuts are applied on the invariant mass of the top 
and Higgs candidate. For both top candidates, the mass 
is required to be in the range of 140 GeV < rrijjj < 



215 GeV. The Higgs candidate mass is required to be in 
the range of 80 GeV < mjj < 150 GeV. 



F. Results 

The estimated signal yields are summarized in 
Tab. IIVI for the case of polarized beams (P e - , P e + ) = 
(—0.8, +0.3), assuming an integrated luminosity of 
1 ab -1 . The resulting distributions for the thrust, Y%^7, 
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TABLE III. Summary of cuts in the analysis of the 6-jet + lepton mode, denoted as 6j. We denote the 4-jet + 2-lepton mode as 
4j, and the 8-jet mode as 8j. Estimated yields are given assuming an integrated luminosity of 1 ab with beam polarizations 
(P e - , P e +) = (—0.8, +0.3). Refer to the text for the details of the b-tagging requirement and the mass cuts. 



ttH(6j) ttH(8j) ttH(4j) tbW ttZ tig* {bb) 

no cuts 2821 28915 68l 980738.5 2406.9 1159.6 

single isolated lepton 179.6 20.7 28.3 340069.0 790.6 397.7 

thrust < 0.77 145.7 18.5 19.2 144999.0 616.7 266.0 

Y s ^4 > 0.005 125.5 16.6 9.2 12297.7 416.2 113.7 

6-tagging 49.0 1.3 2.9 172.9 53.3 37.8 

mass cuts 39.5 1.2 0.4 23.0 33.9 13.2 



tig -^ftbb 
H ttZ 
tbW 

■ ftH(8j/4j) - 
\JtiH(6j) 




120 140 160 180 200 220 240 40 60 80 100 120 140 160 180 
Top candidate mass (GeV) Higgs candidate mass (GeV) 



FIG. 5. The distributions of the cut variables in the 6-jet + lepton analysis are shown: (a) thrust, (b) Ys^4, (c) mass of the 
top candidate, (d) mass of the Higgs candidate. Each sample is weighted assuming an integrated luminosity of 1 ab~ with 
beam polarizations (P e - , P e +) = (—0.8, +0.3). In each of these 4 plots, all the event selection criteria are applied except for 
the cut on the variable shown. The arrows indicate the optimized cut values. 



the top mass, and the Higgs mass are shown in Fig. |B1 
The signal significance in the 8-jet mode is 3.7, corre- 
sponding to the measurement accuracy of the top Yukawa 
coupling of \Ag t /gt\ = 14%. With unpolarized beams 
(P e - , P e +) = (0.0, 0.0), the significance becomes 2.8, cor- 
responding to \Ag t /gt\ = 18%. 

VII. CONCLUSIONS 

We have evaluated the accuracy of the top Yukawa cou- 
pling at = 500 GeV, taking into account the tt bound- 
state effects for the e + e~ — s- ttH signal sample as well as 
the e + e~ — > ttZ background sample. Other backgrounds 
considered were e + e~ — > tbW~ /tbW + — > bW + bW~ and 
e+e~ -> tig* -+ bW + ~b~W~bb. A simple cut-and-count 
analysis was performed for the 6-jet + lepton and 8-jet 
signal decay modes. We assume an integrated luminos- 
ity of 1 ab -1 . Because the 6-jet + lepton sample and 
the 8-jet sample are statistically independent, the com- 



bined significance can be computed simply by summing 
the significances of the two modes in quadrature, assum- 
ing Gaussian statistics. 

With polarized beams (P e -,P e +) = (-0.8,0.3), the 
combined significance is 5.2, corresponding to the mea- 
surement accuracy of the top Yukawa coupling of 
\^gt/gt\ = 10%. With unpolarized beams (P e -,P e +) = 
(0.0,0.0), the combined significance becomes 4.0, corre- 
sponding to \Ag t /gt\ = 13%. Note that these numbers 
only take into account the statistical uncertainty. 

The largest background contribution is the e + e~ — > 
ttZ process which survives the event selection procedure 
primarily because of the overlapping of the dijct mass 
for the Z and the Higgs. This can be reduced by im- 
proving the jet energy resolution and the jet clustering 
procedures, which in turn improves the mass resolution 
of the Higgs candidate. The second largest background 
contribution is the e + e~ -+ tbW process. Thus it will 
be critical to be able to model the e + e~ — > tbW cross 
section accurately, particularly in the tails of its kine- 
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TABLE IV. Summary of cuts in the analysis of the 8-jet mode, denoted as 8j. We denote the 6-jet + lepton mode as 6j, and the 
4-jet + 2-lepton mode as 4j. Estimated yields are given assuming an integrated luminosity of 1 ab with beam polarizations 
(P e - , P e +) — (—0.8, +0.3). Refer to the text for the details of the 6-tagging requirement and the mass cuts. 



tiH(8j) ttH{6j) tiH(4j) tbW ttZ_ tig* (bb) 

no cuts 2891) 2821 68^ 980738.5 2406.9 1159.6 

reject isolated leptons 266.3 85.6 6.6 589716.0 1351.4 701.2 

thrust < 0.7 167.7 44.0 2.7 107227.0 818.0 311.5 

Y 8 ^7 > 0.0009 113.8 13.0 0.3 4048.1 349.6 67.1 

6-tagging 66.6 6.8 0.1 442.6 77.6 39.8 

mass cuts 50.1 0.4 0.0 75.6 47.6 14.1 




Top candidate mass (GeV) Higgs candidate mass (GeV) 



FIG. 6. The distributions of the cut variables in the 8-jet analysis are shown: (a) thrust, (b) Ys->7, (c) mass of the top 
candidate, (d) mass of the Higgs candidate. Each sample is weighted assuming an integrated luminosity of 1 ab -1 with beam 
polarizations (P e - , P e +) = (—0.8, +0.3). In each of these 4 plots, all the event selection criteria are applied except for the cut 
on the variable shown. The arrows indicate the optimized cut values. 



matically allowed region. 

Our results indicate that the measurement of the top 
Yukawa coupling is possible down to the 10% level of 
statistical precision at the ILC with s/s = 500 GeV af- 
ter taking into account the it bound-state effects, which 
agrees with previous predictions [HI . It will be critical to 
reduce the systematic effects down to the level compara- 
ble to the statistical uncertainties. We expect the system- 
atic uncertainties coming from the determination of the 
background rates to be the dominant effect. The amount 
of tbW background can be estimated by measuring the 
tbW cross section at yfs = 500 GeV. The it bound-state 
effects must also be verified by measuring the it cross sec- 
tion at its production threshold (yfs w 350 GeV) which 
will be used to estimate the rate of the e + e~ — > iiH sig- 
nal and the e + e~ — > ttZ background. For this, it will 
be necessary to measure the differential cross section of 
e + e~ -+ tt in order to separate the Higgs-exchange con- 
tribution via the i-channcl which itself contains the top 
Yukawa coupling. 



On the theoretical front, it will be desirable to reduce 
the uncertainties in the production cross section com- 
ing from loop corrections, which will be critical for pre- 
cise background estimation. For the e + e~ — > tt pro- 
cess, the electroweak corrections are known at the 1- 
loop level [28| . with further improvements expected in 
the coming years. QCD corrections are already known 
at the 3-loop level [29M33| . For the e + e _ — > iiZ pro- 
cess, the known QCD corrections at the 1-loop level [2H 
include the tt bound-state effects. Since our study also 
incorporates the it bound-state effects, it will be neces- 
sary to calculate the higher order corrections in order 
to properly estimate the theoretical uncertainties in the 
e + e~ — > itZ cross section. 
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